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Large-scale Structure and Dynamics of the Most X-ray 
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ABSTRACT 

We present photometric, spectroscopic and weak lensing analysis of the large-scale 
structure and dynamics of the most X-ray luminous galaxy cluster known, RX J 1347- 
1145, at z = 0.451. We spectroscopically confirmed 47 new members with LDSS3 on 
the Magellan telescope. Together with previously known members, we measure a new 
velocity dispersion of 1163 ± 97 km s~^. The mass inferred from our velocity disper- 
sion is M200 = 1-I6ta27 X lO^^Af©, with r2oo = 1-85 Mpc, under the assumption 
of a singular isothermal sphere. We also present a weak lensing analysis using deep 
CFHT data on this cluster, and find a deprojected mass of 1.471q'43 ^ IO^^Mq within 
r2oo, in excellent agreement with our dynamical estimate. Thus, our new dynamical 
mass estimate is consistent with that from weak lensing and X-ray studies in the lit- 
erature, resolving a previously claimed discrepancy. We photometrically detect and 
spectroscopically confirm another massive cluster with a — 780 ± 100 km and 
M200 = 3.4tJ;* X IO^^Mq ~ 7 Mpc south-west of RX J1347-1145, which we refer to 
as RXJ1347-SW. Our spectroscopic survey reveals a possible excess of galaxies in ve- 
locity space in the region between RX J1347-1145 and RXJ1347-SW; comparing with 
simulations, this excess appears consistent with that expected from a large filamentary 
structure traced by galaxies connecting these two clusters. 

Key words: Galaxies: Clusters: General, Galaxies: Clusters: Individual (RX J1347- 
1145) 



1 INTRODUCTION 

How galaxy properties are affected by the environment 
they reside in remains one of the unsolved puzzles of 
gala:xy formation. The galaxy population in clusters has 
been observed to have different p roperties from that in 



the field (e.g. iBalogh et al.l |2004| : IWeinmann et al.l I2OO6I : 



iMcGee et al. I2OO8I ). suggesting that this growth of large 



scale structure is accompanied by a transformation in galaxy 
properties. However, it is still an open question as to ex- 
actl y what mechanisms are responsible for this difference 
(e.g. iBalogh et aLlll999l : iTanaka et al]|2006D . 

Clues might come from the outskirt regions of clus- 
ters, because that is where field and groups of galaxies are 
accreting onto clusters, especially along filamentary struc- 
tures that are seen in simulations (e.g. iBond et al. 1 1 19961 : 
IColberg et al ]|l999l '). Groups in the outskirts are particularly 



interesting; a recent study bv lMcGee et al.l |200^, using a 
semi-analytic model, shows that a large fracti on of galaxies 
that accrete onto clusters are in groups (but see lBerrier et al.l 
l2009h . and they might have been pre-processe d in the group 
environment befor e they fall in t o clusters jBa logh et ajj 
I2OOO; Fuiita 2004: iLi et al.l l2009l : IBalogh fc McG cc 200^). 
However, these groups are likely to be different from iso- 
lated groups in the field , due to their earlier formation time 
dMaulbetsch et cil]|2007l ) and the tidal fields fro m the large- 
scale structure around clusters l|Gnedin|[2003l ). Therefore, 
closer examination of groups and filamentary structures in 
the outskirts of clusters, and a comparison with isolated 
groups in the field, may provide important clues to how the 
cluster environment affects the properties of galaxies. 

We have started an extensive spectroscopic campaign 
on two contrasting clusters at intermediate redshift, focus- 
ing on groups in the outskirt regions out to ~ 8Mpc. The 
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ultimate goal is to put strong constraints on the eflFect of 
the environments by comparing properties of groups in the 
infall regions with those of isolated groups, and with galaxy 
formation models. In this work, as the first paper of a series, 
we present results from a pilot study on the dynamics and 
large-scale structure of one of the clusters, RX J1347-1145 
(2=0.451). 

RX J 1347- 1145 is the most X-r ay luminous cluster i n 
the ROSAT All-Sky Survey (RASS, ISchindler et~al]ll995h . 



scopic (|Cohen & Kneib 


2002f). X-rav JSchindler et al.l 19971: 


Allen et al.| 20021, 20081: 


Ettori et al.1 2004 Gitti & Schindlei 


2004: Gitti et al. 200 


p; Bradac et al.l l2008l: Ota et al. 


200g 


: Miranda et al.l 


2008). lensiuE jFischer & Tvson 


1997 


: ISahu et al. 1998 


: iKline et al.1 2005: Bradac et al. 


2005 


. 20081: Halkola et al.1 l2008l: iMiranda et al.1 l2008l: 


Broadhurst et al.' 20081). and Sunyaev-Zel'dovich (SZ) 



effect ( Pointcco uteau et al.l l200ll : iKomatsu et al.l I2OOII : 
iKitavama et al.l |2004| ) analyses. The centre of the cluster 
was initially thought to be dynamically relaxed. Although 
there are two very bright galaxies near the centre (usually 
referred to as the brightest cluster galaxy (BCG) and the 
second BCG), they are only 18 arcsec (~ O.lMpc) apart in 
the plane of sky, and < 100 km/s in velocity. However, a 
discrepancy between the relatively low mass measured dy- 
namically bv lCohen fc KneibI (|2002l ). compared with X-ray 
and weak lensing measurements, led to the suggestion that 
the cluster is actually undergoing a major merger. This hy- 
pothesis has received support recently fro m the discovery of 



shocked gas in th e centre of the cluster ([Allen et al 



Gitti fc Schindlejr2 004: Ot a et al.ll20oi: iBradac et al 



Miranda et al.ll2008i : IKomatsu et al.ll2001 



20021: 



20081 : 



Most of the previous work is on the central region of 
the cluster; in this paper, we focus on the dynamics and 
the large-scale structure of this cluster, using imaging data 
obtained with CFHT and spectroscopic data obtained with 
the Magellan telescope. We describe the data in 321 In 
we present the photometric and spectroscopic analysis on 
the large-scale structure, and the detection of another mas- 
sive cluster, RXJ1347-SW, ~ 7 Mpc from the cluster centre. 
We also present mass estimates of the two clusters through 
dynamical methods and weak lensing analysis. In !j4]we com- 
pare our mass estimates of the main cluster with those in 
the literature, and discuss the possible filamentary connec- 
tion between the main cluster and RXJ1347-SW using sim- 
ulations. Our findings are summarized in iJSl 

Throughout this paper, we assume a cosmology with 
Q.m = 0.3, Q.A = 0.7 and 7^0=70 km s"^ Mpc~\ All magni- 
tudes are in the AB system. 



2 DATA 
2.1 Imaging 

Deep g' , r' ,i' ,z' band images of the cluster over a 1 square 
degree field were obtained using Megaprime on the CFHT in 
semester 05A (proposals 05AC10 and 05AC12), with total 
exposure times of 4200 s in t;', 7200 s in r', 1600 s in i', and 
3150 s in z' . The z' band data will not be mentioned further, 
because they are not relevant for the work presented here. In 
each filter, images from single exposures are co-added using 



SWARF V2.16.0 l|Bertin et al.ll2002l) and photometry is per - 
formed using SExtractor v2.4.4 (jBertin fc ArnoutsI 119961 ) ■ 
with the camera run zero point. The latter quantity is pro- 
vided by Elixir, and is extracted from the Landoldt field 
observed under photometric conditions every night. Since 
our science data were not taken under photometric condi- 
tions, a zero point offset calculated from a short exposure 
taken under photometric conditions is applied. The typical 
seeing is ~ 0.8", ~ 0.7" and ~ 1" in g' , r' and i! band re- 
spectively. Colours are measured using magnitudes within a 
fixed aperture of diameter 4.7 arcsec, which is much larger 
than the seeing difference in different filters. The magnitude 
measured within the Kron radius (mag^auto) is considered 
as the total magnitude for each object. Although crowding 
and blending of objects afi'ects the photometry for about ten 
percent of the galaxies in our sample, we have checked that 
none of the results presented here are sensitive to this. 

Star-galaxy separation is based on the stellar locus de- 
fined by the half-light radius as a function of total mag- 
nitude, in combination with the SExtractor parameter, 
clasa^star. Objects with class^atar < 0.98 in the g' band, 
and not on the stellar locus, are considered galaxies. The 
completeness limit for galaxi es is ~0.8 mag brighter than 
the 5cr limit for point sources (lYeelll99ll) : therefore, the com- 
pleteness limit for galaxies in our catalog in the three bands 
that are relevant in this work is estimated to be g' ~ 25 
mag, r' ~ 25 mag, and i' ~ 23.2 mag. 



2.2 Spectroscopy 

2.2.1 Observations 

In this section, we first describe our spectroscopic observa- 
tion at the four pointings (shown in Figure [l| that are cho- 
sen based on the large-scale structure traced by red-sequence 
galaxies, which will be described in H3.ll 

To make the survey more efficient, we use the combina- 
tion of {g' — r') and (r' — i') colour information to eliminate 
obvious higher redshift galaxies (see figure 6 in ILu et all 
I2OO9I for more details). Moreover, we select galaxies with 
i' < 22 mag. This ensures that for objects of interest here 
((<;' — i') < 3) our sample is complete in each band, and the 
colour measurement is reliable. We make sure our targets are 
not colour biased; as can be seen in Appendix [X] the com- 
pleteness in each pointing for the red and blue populations 
separately is similar, at ~ 50 per cent. 

Spectra were obtained using the Low Dispersion Sur- 
vey Spectrograph 3 (LDSS3) on the 6.5 meter Clay Tele- 
scope over the nights from 2008 April 28th to May 1st. The 
pointing positions are indicated in Figure [TJ each with a 
FOV of 8.3 arcmin in diameter (corresponding to r=1.4 Mpc 
at the redshift of the cluster). The VR filter and Medium 
Red grism were used, which covers the wavelength range 
5OOOA to 72OOA. At the redshift of the cluster, the impor- 
tant spectral features CaH, CaK, G-band, H5 and [OH] fall 
into this wavelength range. Each pointing was observed with 
two masks, with an average of ~ 60 objects and 3x20 mins 
exposure, per mask, with slit width 1.2 arcsec, providing a 
spectral resolution of ~ 14A. 
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Figure 1. The spatial distribution of tlie large-scale structure with the surface density contours of red galaxies overlaid. The density 
contours start at 2cr from the mean and increase with an interval of la. The small circles indicate all galaxies with red colours, the 
squares and crosses are respectively confirmed red and blue galaxies that are associated with the large-scale structure. The four large 
circles indicate the pointing positions with spectroscopic follow-up. 



2.2.2 Data Reduction 

The 2d raw spectra were first reduced using tiie COSMOS2 
software package, then Id spectra were extracted and red- 
shifts were determined using a m odified version of z spec 
kindly provided by Renbin Yan (jCooper et al ] I2OO9I ). In 
zspec, spectra are first compared with eigen-templates and 
the ten best-fitted redshifts are provided. We then visually 



examine the spectra and determine the best redshift. Usu- 
ally the first or second choice provided is a good fit. Occa- 
sionally when none of the ten choices is a good match, we 
manually cross-correlate the spectral features and find the 
best-fit redshift. We consider redshifts determined from at 
least three supporting spectral features as secure. Two ex- 
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ample spectra, near the magnitude limit of our survey, are 
shown in Appendix [Cl 

For the 473 objects targeted, we obtained 400 secure 
redshifts, which are used in the analysis here. Further de- 
tails on the success rate, and the efficiency of our survey 
are presented in Appendix |X] From the objects that have 
duplicate observations, the uncertainty on the redshift is es- 
timated to be ~ 113 km s~^ in the rest frame (see Appendix 
[B|. We also compare the reds hift of the commo i i obje cts in 
our catalogue with that from I Cohen fc KneibI l|2002l '). and 
find no systematic difference between these measurements 
f Appendix IB)) : therefore, in the following an alysis, we add 
the extra 34 redshifts from ICohen fc KneibI l(2002l to our 
sample when applicable. 



3 RESULTS 

3.1 Photometrically Identified Large-Scale 
Structure 

In this section, we describe the large-scale structure traced 
by the red-sequence galaxies. Again, only galaxies with 
i' < 22 mag are used, for reasons discussed in i]2.2.1l At 
the redshift of this cluster {z — 0.45), the (<?' — i') colour 
brackets the 4000A break; therefore we use this colour to 
identify red-sequence galaxies around this cluster, using a 
method ba sed on the Red-sequence Cluster Sur vey tech- 
nique (e.g. iGladders fc Ye3l2000l : ILu et al.ll2009D . We de- 
fine the red-sequence at the redshift of the c l uster using 
known cluster members from ICohen fc KneibI l|2002h that 
have {g' — i') colour measurement in our catalogue. Given 
the small slope and the scatter of the colour, we choose to 
use a red-sequence that has zero slope with a half width 
of 0.15 mag. Assuming the probability distribution of the 
colour of each galaxy is Gaussian, with a being the error on 
the colour measured from SExtractor, we can calculate the 
probabilit y of each galaxy belong ing to the colour slice (for 
details see iGladders fc Yeell200Gl ). The fact that the colour 
error from SExtractor is likely to be underestimated is not 
important here, because the width of the colour slice used 
is much wider than the colour uncertainties. We consider 
galaxies with at least 10 per cent probability as red-sequence 
galaxies, and they are plotted as small circles in Figure [T] 
Our results here are insensitive to the details of how the 
red-sequence galaxies are selected. 

Using the background number counts measured 
from the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) Wide field 1, we can calculate the overdensity 
of red-sequence galaxies. We define the significance of the 
overdensity as the background-subtracted number of red- 
sequence galaxies divided by the r.m.s. of the background 
distribution. The contours of the significance map, smoothed 
over a region of radius ~ 0.5 Mpc with a top hat function, 
are plotted in Figure [TJ starting from 2a, with an interval 
of la. 

The most significant peak (^ 25ct) is the main clus- 
ter (the pointing "Centre"). The number of red-sequence 
galaxies (background subtracted) within 0.5 Mpc from the 
cluster centre, brighter than m* +2 {i' — 22 mag), denoted 
Nred,m*+2, Is ~ 46 ± 7, where the uncertainty is estimated 
using Poisson statistics. Converting Nred,m*+2 to A''2oo ac- 




5 



Figure 2. The galaxy redshift distribution from the full area 
is shown as the dotted, open histogram. The shaded histogram 
shows the distribution of galaxies in the centre of the main cluster 
(r <^ 1.4 Mpc), with a Gaussian function that is statistically con- 
sistent with this distribution, overlaid. The redshift distribution 
of galaxies within r ^ 1.4 Mpc around RXJ1347-SW is shown as 
the solid open histogram. 



cording to figure 8 in iLu et al.l l|2009D and using the scal- 
ing relation between A'200 and M200 (| Johnston et al.ll2007l : 
Hansen et al. 2007), we estimate the mass of the main clus- 
ter from its richness to be M200 ~ l-4lo3 ^ 10^^ M©; al- 
though we note the conversion from Nred,m*+2 to for 
such a high richness system is extrapolated from lower rich- 
ness systems. 

There is another prominent peak with a significance of 
~ 10a, about 7 Mpc (projected) away from the main cluster 
(pointing "SW"). The Nred,m*+2 for this system is ~ 22 ± 
5, which corresponds to M200 = (3.3li g) x IO^Mq. This 
suggests this system is a cluster in its own right, and we will 
refer it as RXJ1347-SW. 

The significance contours also reveal several clumps of 
overdensities in the infall region around the main cluster, 
and an elongated large-scale structure running from NE to 
SW, along a line connecting the main cluster and RXJ1347- 
SW. Four positions along this line are targeted for our spec- 
troscopic follow-up with the Magellan telescope as described 
in H2.2.1I (Note the significant overdensity east of the point- 
ing "Mid" has been targeted with our on-going Gemini pro- 
gram.) 



3.2 Dynamics 

In this section, we present the dynamical analysis of the 
above described photometrically identified structure, based 
on the redshifts obtained. 
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3.2.1 Redshift Distribution and Velocity Dispersion 

In Figure (2] the full redshift distribution of all galaxies is 
shown as the dotted open histogram. There are two main 
peaks in the full distribution: one corresponds to the centre 
of the main cluster, and the other corresponds to RXJ1347- 
SW. We analyze these two clusters in detail first in this sec- 
tion, and examine the dynamics of the large-scale structure 
in the next section. 

The redshift distribution of galaxies within ~ 1.4 Mpc 
from the main cluster centre is shown as the shaded his- 
togram. The distribution is smooth in general; there is an 
indication of two separate peaks, but the distribution is still 
consistent with being drawn from the single Gaussian over- 
plotted on the shaded histogram in Figure (2] In particular, 
a test shows that a subsample drawn from this single 
Gaussian distribution would have a value as large as we 
observed about 80 per cent of the time. The mean redshift 
we measure is 0. 4513 ± 0.0006, consiste nt with the previous 
measurement bv lCohen fc KneibI (|2002l ) . The disper s ion w e 
measure, using the bi- weight method of iBeers et al.l (|l990l ). 
with 3(T clipping applied, is 0.0056 ± 0.0005, which corre- 
sponds to a rest-frame velocity dispersion of 1163±97 km 
s~^. The erro rs arc estimated using Jackknife resampling 
IWall fc Jenkin s 2003). From now on, we refer to galaxies 
that are within 3 times the dispersion as members of the 
main cluste r. The velocity d i spersi on we measure is higher 
than that of lCohen fc KneibI (|2002l ') (we discuss this more m 

ami. 

Despite the consistency with a single Gaussian there is 
an indication of two separate velocity peaks in the distribu- 
tion. A test shows that, at the ~30 per cent level, the dis- 
tribution is consistent with a combination of two Gaussian 
functions: one with zi = 0.4495 and rest-frame velocity dis- 
persion (Ti — 620 km/s, and the other with Z2 — 0.4552 and 
o"2 = 535 km/s. In that case, the centre of the cluster would 
be composed of two separate components with much smaller 
masses, and a relative velocity of ~ 1000 km s~^; consider- 
ably larger than the < lOOkm/s difference between the two 
BCGs {zi = 0.4515 and 22 = 0.4511). Although statisti- 
cally this model is not favoured over the single-population 
model, we cannot rule out the possibility that there are two 
distinct dynamical components. We revisit this further in 
§ 14.1.21 when we consider both the spatial and velocity dis- 
tribution of galaxies. 

The redshift distribution of RXJ1347-SW (pointing 
"SW" ) is shown as the open solid histogram in Figure (2] 
It is not significantly different from a Gaussian distribution, 
with a mean redshift of 0.4708 ± 0.0006, and dispersion of 
0.0038 ± 0.0005, or 780 ± 100 km s'V 



3.2.2 Phase Space Distribution 

We examine the dynamics of the large-scale structure by 
converting the redshift of each galaxy to a rest-frame ve- 
locity relative to the central redshift of the main cluster. In 
Figure [3] we show the velocity offset as a function of po- 
sition projected along the line connecting the main cluster 
and RXJ1347-SW. Galaxies in the inner region of the cluster 
mainly occupy regions spanning from —2000 km s~^ to 2000 
km s"\ symmetrically, while RXJ1347-SW is at ~ 4000 km 
s^^ from the main cluster. Note that in velocity space, the 
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Figure 3. The rest frame velocity offset as a function of distance 
relative to the cluster centre projected along the line towards 
RXJ1347-SW. The dashed curves indicate the boundary of the 
caustics. See text for details. 



region occupied by members of massive clusters resembles a 
trumpet shape, as expected in a spherical infall model (e.g. 
Ivan Haarlem fc van de WevgaertlliggS ). The dashed curves 
in Figure [3] indicate the boundary of the caustics for these 
two systems (as will be discussed in i )3.4.1 1, derived ana- 
lytically from equation (7) in Diafcrio ( 199£ ) assuming the 
two systems follow a NFW profile (|Navarro et al.|[l997l ). In 
pointing "East" , on the other side of the main cluster away 
from RXJ1347-SW, there are almost no galaxies in the re- 
gions labelled "A" and "B" in Figure |3] around the clus- 
ter above the caustic curve. In comparison, there are ~ 8 
galaxies occupying the region between the main cluster and 
RXJ1347-SW, labelled "C". Given that our redshift com- 
pleteness is about 50 per cent, we expect there are approx- 
imately twice as many i' < 22 galaxies (~ 16) within that 
region in total. Note this is not a sampling bias, because 
the pointing "East" and "Mid" are located symmetrically 
on two sides of the main cluster. Also, the uncertainties on 
the redshift determination, or the position of the caustic 
curve, would not change the fact that there is an excess of 
galaxies between the main cluster and RXJ1347-SW. In M.2\ 
we will use simulations to investigate whether these galaxies 
are physically residing in regions between the main clus- 
ter and RXJ1347-SW, and thus suggest a possible structure 
connecting the two systems. 



3.3 Colour Distribution 



As described in gSTT] we trace the large-scale structure 
around this cluster using red-sequence galaxies. We consider 
galaxies in the redshift range Zd — Sad < z < Zsu> + 3o"sm 
as associated with the structure, where Zd and are the 
mean and dispersion of the redshift distribution of the main 
cluster, and Zsm and asw are the same for RXJ1347-SW. In 
Figure [TJ solid squares and crosses represent red and blue 
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Figure 4. The colour-magnitude diagram of galaxies is shown 
for each of the four pointings we targeted, along the elongated 
structure. Open circles are all galaxies in the field of view, crosses 
are galaxies with secure redshifts, and filled circles are members 
associated with the structure. Solid lines are fitted red-sequences, 
and dashed lines separate blue galaxies from the red-sequence (see 
text for detail). 



galaxies that are spectroscopically confirmed to be associ- 
ated with the struct ure, respectively. 

The study by iButcher fc Oemled l|l984h revealed a 
higher fraction of blue galaxies in high redshift clusters com- 
pared to local ones (the Butcher-Oemler effect), indicating 
on-going star formation in high redshift clusters. Since then 
there have been many studies on the dependence of this ef- 
fect on the environment, and cluster propertie s such as the 
X-ray temperature (e.g. lUrguhart et al. I l2009l : IWake et al.l 
I2OO5I ). RX J1347-1145 is the most X-ray luminous cluster 
known, and has an elongated large-scale structure; there- 
fore it is interesting to explore the colour distribution of the 
galaxies in this structure. The colour-magnitude diagram of 
galaxies in the four pointings along the structure is shown 
in Figure |4l Open circles are all galaxies in the field of view, 
crosses are galaxies with secure redshift obtained, and filled 
circles are members associated with the structure. 

W e calculate the blue fracti on of the main cluster us- 
ing the lButcher fc Oemletl (|l984| ) criterion, i.e. galaxies 0.2 
mag bluer than the CMR in rest-frame (B-V). We transform 
rest-frame A{B — V) = 0.2 into our observed A{g' — i') ns- 
ing the colour differe nce between E and Sab galaxies from 
iFukugita etlo] l|l995h . as indicated by the dashed lines in 
Figure |31 The number of blue cluster members is calcu- 
lated by dividing the number of spectroscopically confirmed 
blue galaxies by the completeness of the spectroscopy for 
the blue population. The number of red cluster members is 
calculated in the same way, and the total number of mem- 
bers is then the sum of that of the blue and red galax- 
ies. The blue fraction of RX J 1347- 1145 within the cen- 
tral 1.4 Mpc (~ 0.8r-2oo), down to i' = 21 (corresponding 
to Mv ~ —20) is 20tg%. The uncertainty is estimated us- 



ing the weighted binomial statistics iGehrelsl (|l986l ). Note 
although RX J1347-1145 is the most X-ray luminous cluster 
in the RASS, its blue fraction in the centr al region is sim- 
ilar t o that of other less luminous clusters (|Ellingson et al.l 
l200ll ). consiste nt with the lack of X-ray luminosity depen- 
dence found bv lWake et al. I ([2005). The blue fraction of the 
cluster RXJ1347-SW (361^2%) is, within the uncertainty, 
consistent with that of the main cluster. 

From the contours in Figure [T] the region between the 
main cluster and RXJ1347-SW shows an overdensity of red- 
sequence galaxies. Our spectroscopy also reveals a large pop- 
ulation of blue galaxies, and it is interesting to note that they 
form two trails running in the same direction as the elonga- 
tion of the large-scale structure. The large blue population 
associated with an overdensity of red galaxies is particularly 
interesting. It suggests that targeting regions with overden- 
sities of red-sequence galaxies for spectroscopic follow-up is 
a strategy that allows the efficient selection of cluster mem- 
bers, while still tracing interesting structures with active 
star formation. A more thorough analysis of the stellar pop- 
ulations in the large-scale structure will be done when the 
rest of the data are obtained. 



3.4 Mass Estimate 

In this section, we estimate the mass of the main cluster and 
RXJ1347-SW using different methods. We defer a compari- 
son with the literature, and the discussion of the dynamical 
implication, to 8^4.11 

3.4-1 Dynamical Estimates 

For a given density profile and velocity anisotropy, the veloc- 
ity dispersion can be used to infer the virial mass. Under the 
assumption of a singular isothermal densit y profile and a ve- 
locity anisotropy that satisfies a — Vc/\/2 (|Lokas fc MamonI 
1200 j ). the virial mass can be expressed as M„ = 2(j\rvlG. 
Applying this to the measured velocity dispersion ( ^3.2. l|l 
yields a M200 = ^^2al| {hH{z)G\ = Liet^J-J? x IO^^Mq, 
with r2oo = \/2cr/ [10-ff (z)] = 1.85 Mpc. This mass estimate 
agrees well with that inferred from the richness of the clus- 
ter, as described in i]3.1l 

Under the same assumptions, the mass of RXJ1347- 
SW is estimated to be M200 = 'iA^W x 10^"* Mq, with 
^■200 = 1.22 Mpc. This M200 is also in good agreement with 
that inferred from the richness, and confirms that the clus- 
ter RXJ1347-SW we detected photometrically using the red- 
sequence galaxies is indeed a massive cluster. There has not 
been any report on the existence of this cluster in the liter- 
ature. 

Another method to estimate the mass through dy- 
namical analysis is the caustic technique, under the as- 
sumption that clu sters form through hierarchical accretion 
(IDiaferio fc GelletHl997l : lDiaferiolll999l ). The advantage of 
this technique is that it does not depend on the assumption 
of dynamical equilibrium, and it can probe the mass out 
to large radii where the virial theorem does not apply any 
more. 

As described in H3.2.2I (Figure [5} , cluster members oc- 
cupy a region in velocity space with a trumpet shape delin- 
eated by caustics that separate them from the Hubble fiow. 
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Figure 5. The deprojected mass profile for RX J1347-1145. Left panel: the black points represent the mass profile inferred from the 
caustic method. The solid, and dotted curves show the fitted NFW and isothermal profile respectively. The stars and open squares are 
our weak lensing analysis. Right panel: the black points are the same as in the left panel. The other symbols are mass estimates from 
the literature. 



Assuming spherical symmetry, the mass of the cluster can 
be estimated as: 

M{<r) = ^l^^ A\r)dr, (1) 

where A is the amplitude of the caustic. The prescription 
fo r implern e ntiu K the caustic method is described in detail 
m iDiaferid (|l999l . Here we only summarize the key steps. 
Since the location of the caustics is where clusters are sep- 
arated from the Hubble flow, it is indicated by a significant 
drop of the density of cluster members in velocity space. 
Therefore, one first calculates the density as a function of 
velocity at different radii, using an adaptive kernel method 
(Silverman 1986). In order to use the spherical smoothing 
window, the window size in velocity and distance, hv and hr, 
has to be rescaled (see Fpiafcrio 1999, for details). For our 
data, we choose q = hv/h,- = 6, and a spherical smoothing 
window size of ~ 700 km/fl In the central region of clusters, 
the escape velocity (essentially the location of the caustics) 
is twice the velocity dispersion, which is valid under the as- 
sumption that the central region of the cluster is virialized. 
Therefore, in practice, the threshold density where the sig- 
nificant drop occurs (and thus determines the amplitude of 
the caustics) is chosen so that the above condition is satis- 
fied. 

The deprojected mass profile of the main cluster esti- 
mated using this technique is shown as solid dots in both 
panels of Fi gure [5l The errors are estimated using Jackknife 
resampling l|Wall fc Jenkinsll2003l '). At r2oo = 1-85 Mpc the 
enclosed mass is ~ 1 x 10^^ Mq. This is in excellent agree- 
ment with that inferred from the velocity dispersion; how- 

^ Within the range of reasonable choices of parameters, the 
change of the resulting mass profile is within the resulting error 
bars. 



ever, this is largely true by construction, since the density 
threshold is chosen such that the two techniques agree within 
the virialized region of the cluster. The caustic-based mea- 
surements have the most value at larger radii, where the 
assumption of virial equilibrium is unlikely to hold. Note 
beyond r ~ 2 Mpc, we only have sparse sampling in two 
pointings, and the velocity distribution of galaxies in these 
pointings is not likely representative of that of the whole 
population at that radius, due to the presence of elongated 
large-scale structure. We fit the mass profile to a NFW 
(Navarro ct al. 1997) and isothermal profile, shown as the 
solid and dotted curves respectively in the left panel of Fig- 
ure!^. Because of the large uncertainties on the mass profile 
due to limited data, the two profiles cannot be distinguished. 

3.4- S Weak Lensing Analysis 

We use the deep CFHT r' imaging data to study the cluster 
mass distribution using weak gravitational lensing. Unlike 
dynamical methods, the weak lensing signal provides a di- 
rect measurement of the projected mass, irrespective of the 
dynamical state of the cluster. Furthermore, the observed 
sign al can be used to "map" the distribution of dark matter 
(e.g. iKaiser fc Sauireslll993l ). 

Our cluster w eak lensin g analy sis follows in detail the 
steps described in iHoekstral (|2007l l . In order to relate the 
lensing signal to a physical mass requires knowledge of the 
redshift distribution of the sources. The wavelength cover- 
age and limited depth of the observations prevent us from 
deriving accurate photometric redshifts for the ind i vidua l 
sources. Instead we use the results from lllbert et al] (|2006l ) 

^ The fitted NFW profile is used in the analytic derivation of the 
caustic curves shown in Figure [3] and discussed in i|3.2.2l 
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Figure 6. The grey scale shows the weak-lensing signal-to-noisc 
map, out to ~ 10 Mpc from the main cluster centre. Contours 
of the red galaxy surface density (as in Figure [T]l are overlaid, 
starting from 2a with an interval of Itr. The two prominent peaks 
are the main cluster and RXJ1347-SW; their mass concentrations 
coincide with the overdensities of the red-sequence galaxies. 



who determined photometric redshifts for galaxies in the 
CFHTLS deep fields. 

To determine the lensing signal, we measure the shapes 
of galaxies with 22 < r' < 25. We use such a rela- 
tively faint cut because of the high redshift of the clus- 
ter. The galaxy shapes are co rrected using the method pre- 
sented in iKaiser et all (Il995t) . with modifications described 
in iHoekstra et al.l ( 19981 ). The weak lensing signal is sensi- 
ti ve to all s t ructure alo ng the line-of-sight and, as shown 
in iHoekstral (|200ll . [20031 ) . this distant (uncorrelated) large- 
scale structure introduces additional noise to the mass es- 
timate. This 'cosmic shear' noise is included in the errors 
listed below. We note that the noise due to the intrinsic el- 
lipticities of the sources dominates because of the relatively 
high cluster redshift. 

The weak lensing reconstruction of the projected sur- 
face density with a smoothing scale of 123.3" is presented in 
Figure out to ~ 10 Mpc. To indicate more clearly which 
structures are significant, we show the signal-to-noise ratio, 
rather than the signal itself, as the grey scale (see the left 
panel in Figure|9]for a more detailed view). The main cluster 
and RXJ1347-SW are clearly detected, at ~ 8(t and ~ Act 
levels respectively. No other significant structures are evi- 
dent in this map. The contours of the overdensity of the red- 
sequence galaxies, the same as in Figure [TJ are overplotted 
on top of the mass map. For the main cluster and RXJ1347- 
SW, the mass concentrations coincide with the overdensities 
of red-sequence galaxies. 




RXJ1347-SW 



r [arc 



Figure 7. Tangential distortion as a function of the distance from 
the cluster center for the main cluster (top panel) and RXJ1347- 
SW (bottom panel). In both cases we detect a significant lensing 
signal out to large radii (the arrows indicate 2 Mpc). The open 
points correspond to the signal that is measured when the sources 
are rotated by 45 degrees, which provides a measure of the level of 
systematics. The points are consistent with zero. The blue dashed 
line shows the best fit SIS model, where the solid blue part indi- 
cates the region to which the model was fit. The error bars indicate 
the statistical error in the measurement, due to t he intrinsic el- 
lipticities of the sources (see iHoekstra et al.ll200Cll for details) . 



The smoothing of the mass map limits its use to deter- 
mine the cluster mass. Instead, we compute the mean tan- 
gential distortions as a function of distance from the cluster 
center. The resulting measurements for the main cluster and 
RXJ1347-SW are shown in Figure[7l We estimate the cluster 
masses from these shear measurements using two different 
methods. The first, and most commonly used, is to assume 
a density profile and determine the best fit parameters to 
the observed tangential distortion. For comparison with the 
dynamical results, we first fit a singular isothermal sphere 
model (SIS) to the observations between 500 kpc and 2 Mpc, 
which yields a velocity dispersion of ct = 1282±124 km/s, in 
excellent agreement with the results based on the redshifts 
of the cluster galaxies. Thanks to the wide-field lensing data, 
we can also apply the same method to RXJ1347-SW (cen- 
tred at 13''46"27.2", -11°54'20.2"), for which we obtain a 
velocity dispersion of ct = 927 ± 165 km/s, also in good 
agreement with the value inferred from the redshifts of the 
cluster members. 

Although the SIS allows for an easy comparison to the 
observed velocity dispe rsion, a more physi cally motivated 
function is the NFW (jNavarro et al.1 Il997l ) profile, which 
has been inferred from numerical s imulations. We ass ume 
the mass-concentration relation from 'Duffy et alj (|2008l ) and 
fit the NFW model to the data within 500 kpc and 2Mpc. 
With these assumptions, we find a best-fit value for the 



^ The map presented in this Figure was constructed using source 
galaxies with 23 < r' < 25.5 



virial mass of Mv 



X W^Mq. The best fit 



NFW model within different apertures are indicated by stars 
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* Fischer & Tyson 97 (weak lensing) 
■ Gitti et al. 07 (X-ray) 
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(strong lensing) 
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□ Kling et al. 05 

(weak lensing) 



Fp (Mpc) 

Figure 8. The projected mass as a function of radius for RX 
J1347-1145. The solid hexagons are estimates from this work, and 
the other symbols are from the literature. See text for details. 



in the left panel of Figure [S] For RXJ1347-SW we obtain 
Mvir = &Atll X 1O"M0. 

To compare to the results from the cluster galaxy dy- 
namics, we also compute the corresponding masses within 



,.dyn 



r 200 — 1-85 Mpc. For this choice of radius we obtain M200 = 
l-47lo.l3 X 10^^ Mq for the main cluster. This result is in 
excellent agreement with the value inferred from dynamics. 



For RXJ1347-SW we obtain M200 
».dyn 



3 + 2.9 



X W^Mq for 



200 - 1.22 Mpc. 

Since weak lensing signals are sensitive to the total mass 
projected along the line-of-sight; it is interesting to consider 
a non-paramet ric estimate for the projected mass withi n an 
aperture (e.g. IClowe et al. I Il998l : iFahlman et al. 1 11994 see 
Hoekstra (2007) for our implementation). The resulting pro- 
jected masses as a function of radius are presented in Figure[S] 
as solid hexagons. Figure [S] also shows other estimates from 
the literature, which will be discussed in det ail in 8^4.1.11 

As discussed in Hoekstra (2007) and iMahdavi et al.l 

(|2008l ) one can assume spherical symmetry to deproject the 
resulting masses. We assume an NFW profile to extrapolate 
along the line-of-sight, but note that the results at large 
radii do not depend strongly on this assumption. The re- 
sulting deprojected masses are shown as open squares in the 
left panel of Figure (5) They are consistent with our dynam- 
ical mass estimate using the velocity caustics, within the 
uncertainties. 

Within an aperture of rjoo = 1-85 Mpc, we find a 
deprojected mass for the main cluster of M200 = 1.25 ± 
0.19 X lO^'^Mo, again in excellent agreement with the dy- 
namical measurements. Similarly we find a mass of M200 = 
5.6±1.6 X 1O"M0 for RXJ1347-SW using = 1.22 Mpc, 
which is somewhat larger than the dynamical mass. 



4 DISCUSSION 

4.1 Dynamical State of the Main Cluster 

The mass estimate and the dynamical state of this cluster 
have been the focus of many previous studies. In this sec- 
tion, we first compare the mass estimated using different 
methods in this work and in the literature, and then discuss 
the dynamical state of the cluster. 

^.1.1 Comparison of Mass Estimates 

The mass of RX J1347-1145 has been estimated using differ- 
ent techniques. Other than through the dynamics and weak 
lensing signals as described in i]3.4l X-ray data have also 

been used to infer the mass. 

In a recent study by iGitti et al.l l|2007l ) , they estimate 
the m ass by finding a NFW mass profile (jNavarro et al.l 
Il997h that reproduces the observed X-ray temperature pro- 
file. The resulting integrated mass at different radii are plot- 
ted as open triangles in the right panel of Figure (5] This is 
consistent with our estimates from both the dynamics and 
the weak lensing analysis ( ii3.4l) 



Similarly, the study by Allen et al.l (|2002l ) estimated the 
virial mass by fitting a NFW profile to their Chandra ob- 
servations, and their Af200 is plotted as the solid hexagon in 
the right pa nel of Figure [5] Other mass estimates us i ng X - 
ray data by ISchindler et al.l l|l997h and lEttori et aP (|2004l ) 
are plotted in the right panel of Figure [5] as well. All these 
estimates at different radii are consistent within their error 
bars. 

As discussed in tj3.4.2l the projected mass is a more di- 
rect measurement than the deprojected mass from the weak 
lensing analysis, since it is the total projected mass along 
the line-of-sight that gives rise to the weak lensing signals; 
therefore, it is more direct to compare the measured pro- 
jected mass from different studies. 

In Figure [S] the projected mass enclosed within a cer- 
tain radius estimated from strong, weak, and combined 
lensing analyses in the literature are plotted, along with 
our weak lensing estimate (see legends). Our weak lensing 
mass estimate (solid h exagons) is consistent with that of 
iFischer fc TvsonI |l993) (the solid triangle), and these two 
estima tes are them selves in agreement with the X-ray anal- 
ysis of iGitti et al.l |2007() (the solid squa re) . 

The weak lensing mass estimate bv lKUng et al] (|2005l ) 
is slightly higher (~ 2a). This discrepancy is likely due 
to differenc e s in t he assumed source redshift distribution; 
iKling et al.l l|2005h base their estimate on r < 24.5 spec- 
troscopy over the small area of the Hubble Deep Field, 
while we use the photometric redshift distribution from the 
CFHTLS Deep survey. Given our good agreement with both 
X-ray and dynamical mass estimates, it seems likely that 
their analysis underestimates the mean source redshift, and 
thus overestimates the cluster mass. 

In c omparison, the mas s estimates from stron g lensing 
studies (|Halkola et all [2008'; 'Mir anda et"al ] '2008") and the 
combi ned strong and weak lousing analysis by Brada 5~et al.l 
(|2008l ) are systematically higheijj In general strong lensing 



* We did not include the work bv lSahu et al.l lll998l ) because the 
redshift of the arc they assumed is not correct. 
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provides an accurate estimate of the projected mass only 
within the central few hundred kpc, whereas the weak lens- 
ing mass is better constrained at larger radii; therefore cau- 
tion should be taken when combining strong and weak lens- 
ing analysis. See also Halko la et al. ( 2008) for a summary 
and more discussion of the comparison among these works. 

The good agreement shown in Figures [5] and |8] suggest 
that the mass estimated from dynamical, X-ray, and weak 
lensing studies is generally consistent, while strong lensing 
gives a slightly higher value. 

So far, there has bee n only one other publish ed redshift 
survey on this cluster, bv lCohen fc KneibI (|2002D . They de- 
rived from 47 cluster members in their sample a velocity 
dispersion of 910±130 km s~ ^, using the same bi- weight 
method by iBeers et al.l (Il990l ) as used in this work. This 
is lower than the value we derived, and thus leads to a lower 
mass that is in conflict with that from other methods. How- 
ever, we believe that the value quoted in their paper is in 
error. When we reanalyze the 47 redshifts as given in their 
Table 1, we obtain a velocity dispersion of 1098 ± 157 km 
s~^. Therefore, there is actually no discrepancy between the 
dynamical mass estimate and that from X-ray and lensing 
analysis. 



4-1.2 Merger in the centre? 

One important reason the merger scenario was proposed 
for this cluster is to explain the large discrepancy be- 
tween the mass estimated from the velocity dispersion by 
ICohen fc KiieibI (|2002l ) and the X-ray and lensing mass es- 
timates discussed in the above section. However, as shown 
above, our new measurement resolves this discrepancy, and 
thus removes the original motivation for believing the sys- 
tem is undergoing a major merger. 

Nonetheless, there is other evidence suggesting that 
this cluster might have recently undergone a merger. Al- 
though the X-ray morphology is generally relaxed in the 
centre of the cluster, the existence of a region with shocked 
gas ~ 10—20 arcsec (~ 0.09 Mpc) SE of the cluster cen- 
tre has been detecte d as a high-pressu re region in X- 
ray observations ( Allen et al.ll2002l: ICitti fc Schindleill200 ' 



20081 : iMiranda et al.H2008l l 



lota et all l2008l : iBradac et al 
and a s trong decrement in t he Sunyae v-Zel'dovich obser- 
vations (jKomatsu et aLlbOOll ). Also, Mir anda et"all (120081 ) 
found that their strong lensing data can be much better de- 
scribed by a two component model with comparable mass 
than the single mass component model. The fact that there 
are two very bright galaxies in the centre of this cluster 18" 
(~ 0.1 Mpc) apart in the plane of sky is also taken as a 
possible indication of a merger event (but n ote their line-of- 
sight relative velocity is only ;Sl00 km/s). iKitavama et al.l 
()2004l ) suggested that one possible scenario is a collision of 
equal mass clumps with a collision speed of 4600 km s~^. 
However, this collision cannot be mostly along the line-of- 
sight, because as sho wn in 8^3.2.11 (and also in the work by 
ICohen fc Kneibll2002l ). the redshift distribution is better de- 
scribed as a single Gaussian distribution. As discussed in 
8^3.2.11 our velocity data are also consistent (at the 30% sig- 
nificance level) , with a model distribution having two Gaus- 
sian components. However, even in that case, the radial ve- 
locity offset between the two peaks is only about 1000 km 
s~^, much smaller than 4600 km s^^. Furthermore, if the 



centre of the cluster is indeed composed of two merging com- 
ponents, each with a much smaller velocity dispersion than 
what we find for the system as a whole (see i l3. 2.1(1 . the 
dynamically-estimated mass profile would be much lower 
than that shown in Figure [5] (the hexagon and the filled 
points), and that would be inconsistent with the weak lens- 
ing and X-ray mass estimates. 

To investigate this further, we use the A-test proposed 
by iDressler fc ShectmanI (|l988l ) to examine closely the 3D 
substructures in the centre of the main cluster. This test 
measures the deviation of the local velocity and velocity dis- 
persion from the global mean, and is quantified as : 



^2 ^ iVnft + 1 



[{Vlocal - + {Olocal - o)^] 



(2) 



where N^b is the number of neighbour galaxies, Viocai and 
(yiocai are local mean velocity and dispersion measured from 
each galaxy and its A'„i, nearest neighbours, and v and cr are 
global mean velocity and velocity dispersion, respectively. 

In the left panel of Figure |9] we show the spatial distri- 
bution of the confirmed cluster members (solid squares) with 
their local deviation from the global mean indicated by the 
size (proportional to 5) of the circle around them. The solid 
circles indicate negative local deviations {z < 0.4513), and 
dashed open circles positive local deviations {z > 0.4513). 
To examine closely the centre of the cluster, in this panel 
we only plot the region within a radius of ~ 0.5 Mpc from 
the cluster centre. The two triangles indicate the two BCGs. 
The region ~ 20 arcsec SE of cluster centre, indicated by the 
large green filled circle, is where the shocked gas is detected. 
Around that region, we do not detect significant deviation 
from the global mean. However, there seems to be a hint of 
a segregation between the region NE and SW of the clus- 
ter centre, in the sense that in the NE region most galaxies 
have negative local mean velocity relative to the mean red- 
shift of the cluster, while galaxies in the SW region mostly 
have positive local mean relative velocity. This segregation 
is along the same direction connecting the main cluster and 
RXJ1347-SW. 

Note the above estimated local mean might be biased 
by a few galaxies with large deviations due to our sparse 
data sampling. Therefore, in the right panel of Figure we 
examine directly the spatial distribution of the redshift of 
individual galaxies, with a field of view of ~2 Mpc x 2 Mpc. 
The open circles indicate possible potential members (not 
targeted), and the crosses represent nonmembers. Cluster 
members with z < 0.4513 are indicated by the blue solid 
squares, and members with z > 0.4513 are indicated by 
the red solid squares. There is a hint that more galaxies 
with negative relative velocity reside in the region NE to the 
cluster centre. To test whether the distribution of the two 
samples is significantly difi'erent, we apply a 2-dimensional 
Kolmogorov-Smirnov test following the procedure proposed 
by Peakcock (1983) . The test suggests that there is a ~60 
per cent probability that galaxies with negative and positive 
relative velocity have different spatial distributions. 

The above evidence discussed in this section suggests 
that this might be a merger more in the plane of the sky 
than along the line-of-sight. However, with the lack of both 
a clear sharp shock front and a distinct separation of mass 
concentration on both sides of the shock (see contours in the 
left panel of Figure [SJ , as observed in the famous example 
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Figure 9. Left panel: the spatial distribution of galaxies in the centre of the main cluster (the field of view is ~1 Mpc X 1 Mpc), with 
the contours of the weak lensing signal-to-noise mass map overlaid, with an interval of 0.5a-. For clarity, we only plot the confirmed 
cluster members (the solid squares). The two triangles indicate the two BCGs. The size of the open circles (proportional to 5) around the 
cluster members show the local deviation from the global mean, with solid open circles indicating negative local deviations (z < 0.4513) 
and dashed open circles positive local deviations (z > 0.4513). The big filled green circle indicates the region ~ 20 arcsec SE from 
cluster centre, where the shocked gas is detected. See text for details. Right panel: the spatial distribution of galaxies in the centre of 
the cluster, with a field of view of ~2 Mpc X 2 Mpc. The open circles indicate possible potential members (not targeted), the crosses 
indicate nonmembers, the blue solid squares represent confirmed members with 2 < 0.4513, and the red solid squares represent confirmed 
members with z > 0.4513. Again the two triangles indicate the two BCGs. 



of merger in the plane of the sky, the B ullet cluster (e.g. 
iMarkevitch et al.ll2002l : IClowe et allbood ) , the exact merger 
scenario still needs closer examination both in observations 
and simulations. 

4.2 A possible filamentary connection 

As seen in Figure |31 cluster RXJ1347-SW is ~ 7 Mpc (pro- 
jected in the plane of sky) away from the main cluster, and 
has a rest-frame velocity offset of ~4000 km s~^ relative to 
the main cluster. This could be a result of two scenarios. 
One is that RXJ1347-SW is at nearly the same distance as 
the main cluster, and is falling into the cluster at a speed 
of ~ 4000 km s"^ The other is that it is ~ 25 Mpc in the 
background moving with Hubble expansion. We resort to 
simulations to try distinguishing these two scenarios. Since 
what we are looking for is how likely it is for two clusters Uke 
RX J1347-1145 and RXJ1347-SW to have a relative pecu- 
liar velocity as large as observed, we use the Hubble volume 
simulations (Evrard et al. 2002) because of its large box size 
(3000^/i-^Mpc^). 

From the cluster catalogue in the full-sky light cone with 
ACDM cosmology (O = 0.3, A = 0.7, erg = 0.9), we search, in 
a redshift range close to the redshift of RXJ1347, for a pair of 
clusters that have similar masses (Ml ;> 1 X 1O^''M0,M2 ;> 
1 X 10^'* M0), projected distance (~ 10 Mpc), and rest-frame 
velocity offset as RX J1347-1145 and RXJ1347-SW. We do 
not find any system where they are at the same distance, but 



with an observed velocity offset of ~ 4000 km s~^. This does 
not mean definitely that no system could have such a large 
relative peculiar velocity, because this is based on only one 
simulation with a large but still limited volume, under the 
assumption of one particular set of cosmological parameters. 
Among these parameters, erg is the most uncertain one, and 
in fact a lower value than used in the simulation is preferre d 
by the latest WMAP measurement l|Dunklev et al.ll2009l ). 
However, in that case, the peculiar velocity between clus- 
ters would be even lower due to the lower attracting force 
from the lower mass density. Even though we cannot reach a 
definitive conclusion, this does suggest that the probability 
that RXJ1347-SW is falling into the main cluster with such 
a high velocity is low; instead, it is more likely to be in the 
background. 

However, just because RXJ1347-SW is far in the back- 
ground does not necessarily mean that there is no connection 
between it and the main cluster. The fact that we observed 
an excess of galaxies between the two in velocity space (Fig- 
ure [5]) suggests the re might be a physical connection (see 
iTanaka et al]|2009l for a report of the spectroscopic detec- 
tion of a large-scale galax y filament associate d with a mas- 
sive cluster at z — 0.5; also lEbeling et al]|2004l ). We compare 
our observation with simulations to see if the observed ex- 
cess of galaxies could represent a real large-scale filamentary 
structure connecting the main cluster and RXJ1347-SW, 
rather than foreground and background galaxies projected 
in velocity space. To examine this on the scale of individ- 
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ual galaxies, we make use of the galax y catalogue from a 
semi-analytic model by lFont et all (l2008lrl based on the high 
resolution Millennium simulation (jSpringel et al.l[2005l ). Be- 
cause of the small volume of the simulation, there are only 
four clusters that have mass greater than 1 x 10^^ M©. We 
search around these four clusters for a second cluster with 
lower mass and at a similar distance projected in the plane 
of sky and along the line-of-sight from the more massive 
cluster as in the case of our observed RX J 1347- 1145 and 
RXJ1347-SW. For the four pairs of systems we found, we 
calculate the line-of-sight velocity offset of each galaxy rel- 
ative to the more massive cluster, and plot it as a function 
of distance, projected along the direction towards the less 
massive one. For direct comparison with our observed data, 
we only plot galaxies in a strip with a width similar to the 
radius of the FOV as our spectroscopic observation, down to 
the same magnitude limit as our data. In Figure[lO]we show 
one extreme case, where the two clusters (Mi ~ 3x IO^^A^q, 
M2 ~ 2 X 10^^ Mq) are connected by a physical, large-scale 
filament traced by galaxies, with two other clusters in be- 
tween. This structure is clearly visible in velocity space. The 
triangles indicate galaxies that are spatially within a cylin- 
der with a radius of 3 Mpc (comparable to the virial radius 
of the more massive cluster) connecting the two clusters. 
This is an extreme case where the galaxy filament is very 
impressive due to the presence of two other massive systems 
in between the two clusters. However, even in other cases in 
the simulation where the excess of galaxies between two clus- 
ters is smaller (~ 5— 10) than this one, the excess is still due 
to galaxies that physically reside in the region between two 
clusters, instead of being purely due to a projection effect. 
Although this is not definitive, it does support the possibil- 
ity that the excess of galaxies we observe in velocity space 
between the two clusters refiects a real physical structure. 



5 SUMMARY 

In this work, we presented photometric and spectroscopic 
analysis on the dynamics and large-scale structure around 
the most X-ray luminous cluster RX J1347-1145, and found 
the following results: 

• We identified an elongated large-scale structure around 
the main cluster, photometrically using red-sequence galax- 
ies, and discovered a rich cluster (RXJ1347-SW) ~ 7 Mpc 
(projected) away from the main cluster. 

• The velocity dispersion of the main cluster is 1163±97 
km s"\ corresponding to M200 = 1.16t°;27 IO^^Mq, with 
r2oo = 1.85 Mpc. 

• Our dynamical mass estimates based on both the velocity 
dispersion and velocity caustics are consistent with our weak 
lensing analysis, and the X-ray studies in the literature. We 
also find excellent agreement between these mass estimates 
and that estimated photometrically from the abundance of 
red-sequence galaxies (the richness). 

• Our new measurement resolves a previous claimed dis- 
crepancy between_th«_dyii^^ mass estimate and other 
methods ICohen fc KneibI lj2002l ). 



^ Note for our purpose here the exact semi-analytic model used 
is unimportant. 
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Figure 10. This figure shows the rest-frame velocity offset of 
galaxies around two simulated clusters, as a function of distance 
projected in the direction connecting the two. The two clusters 
(Ml = 2.7 X W^^Mq, M2 = 1.8 X 10^5 M0) are connected by 
a real physical filament traced by galaxies, with two other clus- 
ters in between. The excess of galaxies between them in velocity 
space is clearly visible. The two squares indicate the centres of 
the two clusters, and the triangles indicate galaxies that are spa- 
tially within a cylinder with a radius of 3 Mpc connecting the two 
clusters. 



• Our spectroscopic analysis confirmed that the photomet- 
rically identified structure RXJ1347-SW is indeed a massive 
cluster, with a velocity dispersion of 780 ± 100 km s"'^. 

• Our weak lensing analysis shows good correspondence 
between the mass concentration and the overdensity of 
red-sequence galaxies in the main cluster and RXJ1347-SW. 

• Comparing with simulations, the observed excess of galax- 
ies between the main cluster and RXJ1347-SW relative 
to other regions around the main cluster in velocity space 
suggests a possible filamentary structure connecting these 
two clusters. This structure has a large population of blue 
galaxies associated with the overdensity of red-sequence 
galaxies, suggesting that targeting overdensities of red- 
sequence galaxies for spectroscopic follow-up allows the 
efficient selection of cluster members, while still tracing 
interesting structures with active star formation. 

This work is a first result from our large on-going spec- 
troscopic campaign. In future work, we will use the technique 
presented in this work to search and examine the properties 
of other structures, especially groups in the infall region, 
with the on-going data collecting. 
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APPENDIX A: REDSHIFT COMPLETENESS, 
SUCCESS RATE AND EFFICIENCY 

In Figure lAll and IA21 we show the sampling completeness 
(number of galaxies targeted / number of all targets), and 
overall completeness, i.e. the success rate (number of galax- 
ies with secure redshift obtained / number targeted) mul- 
tiplied by the sampling completeness, for the red and blue 
population in the 4 pointings, respectively. Our overall com- 
pleteness is about 50 per cent. 

For red galaxies, in pointings "Centre" and "SW", 
about 90 per cent of galaxies with secure redshifts obtained 
turned out to be members; and in the rest of the two less 
densed pontings, the efficiency is ~60 per cent and ~70 per 
cent. For blue galaxies, the typical efficiency is about 30 per 
cent. (The relatively low efficiency of obtaining blue mem- 
bers is not surprising.) 



APPENDIX B: REDSHIFT UNCERTAINTY 

We possess duplicate observations for some of the objects 
that appear on both masks at each pointing. To determine 
the uncertainty on the redshift, we compare the difference 
in the redshifts of the same object from two observations in 
Figure lBl] The histogram is the distribution of the difference 
in redshift (top axis) or rest-frame velocity (bottom axis) 
from repeated observations. The vertical solid line is the 
mean of the distribution, and vertical dashed lines indicate 
la. It shows that the uncertainty on our redshift is ~113 km 
s~^ in the rest frame. 

We also compare the redshifts determ i ned f rom our 
observation with that from ICohen fc KneibI l|2002l ) for the 
matched objects. The distribution of the difference is shown 
in the left panel of Figure IB2I with top axis in unit of red- 
shift and bottom axis rest-frame velocity. In the right panel 
of Figure IB2I the distribution of the difference in terms of 
the rest-frame velocity is normalized by the quadratically 
a dded uncerta inties on our redshift (113 km s~^) and that 
of lCohen fc Kncib (2002) (which they assumed to be 100 km 
s~^). This distribution is not significantly different from a 
standard normal distribution shown as the curve, although 
it is a small number statistics. This suggests that there 
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Figure Al. The completeness of spectroscopic follow-up and suc- 
cess rate for red galaxies in the 4 pointings. The dashed lines 
show the sampling completeness (number of galaxies targeted / 
number of all targets), and the solid lines show the overall com- 
pleteness, i.e. the success rate (number of galaxies with secure 
redshift obtained / number targeted) multiplied by the sampling 
completeness. 
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Figure A2. The same as Figure IXT] but for the blue population. 



is no systematic offset between the redshifts measured by 
ICohen fc KneibI (l2002l) and us; therefo re, we add the extra 
34 redshift fa-om lCohen fc KneibI (|2002f ) to our sample when 
applicable. 
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Figure Bl. The distribution of the difference in redshift (top 
axis) or rest-frame velocity (bottom axis) for galaxies with re- 
peated observations in our sample. The vertical solid line is the 
mean of the distribution, and the vertical dashed lines indicate 
la dispersion of the distribution. 




Figure B2. Left panel: the distribution of the difference in red- 
shift (top axis) and rest-frame velocity (bottom a xis) for common 
object s between our sample and the sample of ICohen fc KneibI 
(|2002l'l . Right panel: the distribution of the difference in the 
rest-frame velocity for the common objects normalized by the 
quadratically added u ncertainties on our redshifts and those of 
ICohen fc KneibI tOOt i. The curve is a standard normal distribu- 
tion. 
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Figure CI. Example spectra (the black line) of two typical galax- 
ies with i'=21.4 mag and j'=21.8 mag, close to the magnitude 
limit of our sample. The cyan line shows the template spectrum. 
The red line shows the sky spectrum, arbitrarily shifted and renor- 
malized, to indicate where bright sky lines give rise to residuals. 
The positions of different spectral features are indicated with ver- 
tical lines, and in both spectra there are multiple absorption fea- 
tures identified, ensuring a secure redshift. 



APPENDIX C: SAMPLE SPECTRA 

Figure [CT] shows example spectra of two typical absorption- 
line galaxies with i'=21.4 mag and i'=21.8 mag, close to the 
magnitude limit of our sample. Multiple absorption features 
are visible, including prominent CaK, CaH, H5, H/3 lines. 
See figure caption for description of difi'erent lines. 
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